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ABSTRACT 

We measure the Petrosian structural properties of 33 brightest cluster galaxies (BCGs) 
at redshifts z ^ 0.1 in X-ray selected clusters with a wide range of X-ray luminosities. 
We find that some BCGs show distinct signatures in their Petrosian profiles, likely 
to be due to cD haloes. We also find that BCGs in high X-ray luminosity clusters 
have shallower surface brightness profiles than those in low X-ray luminosity clusters. 
This suggests that the BCGs in high X-ray luminosity clusters have undergone up 
to twice as many equal-mass mergers in their past as those in low X-ray luminosity 
clusters. This is qualitatively consistent with the predictions of hierarchical structure 
formation. 

Key words: Galaxies: clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: evolution - galaxies: formation 



1 INTRODUCTION 

The first-ranked or brightest cluster galaxies (BCGs) are 
uniquely positioned at the centre of the gravitational po- 
tential of clusters of galaxies. With brightnesses roughly ten 
times that of L* galaxies, BCGs are not simply t he bright 
exten s ion of a Schcch ter luminosity f unction (e.g. ISandage 
; iTremaine fc Richstond Il977l iBernstein fc Bhavsar 
20011) . They are predominantly elliptical in their morphol- 



ogy and owing to their extreme luminosities and position, 
they dominate their host cluster visually, forming a unique 
population in their own right. However, their evolution and 
how their environment affects that evolution, are still poorly 
understood. 

In previous papers of this series on the K-ban d 
photometric prog erties of BCGs; Collins & Mannj il99d) ; 
iBurke et al] feoool) and iBrough et all tOQ-j) : we observe 
that BCGs in clusters with an X-ray luminosity Lx > 
1.9 X 10'''' ergs 

= 0.7, nm 

form than those in their low X-ray lu minosity count e rparts 
at redshifts z ^ 0.1. It appears from iBrough et al.l ||2002D 
that the evolutionary history of BCGs depends on the rich- 
ness of their cluster environment, with BCGs in the most 
luminous clusters showing little sign of merging since 2; ~ 1 
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^ (in the EMSS passband 0.3-3.5 keV with 
- 0.3, h = 0.7) are brighter and more uni- 



and those in the least luminous clusters showing significant 
signs of mass growth over the s ame timescale . Howe ver, the 
mass accretion rates inferred in lBrough et al.l ll2002ll rely on 
the assumption that increases in luminosity are attributable 
to increases in mass. 

It is important to confirm this result to determine 
whether the mass evolution of BCGs is dependent on their 
environment. If BCGs are in dynamical equilibrium then 
changes in mass should produce changes in size. Therefore, 
a means to test the assumption of increases in luminosity be- 
ing attributable to mass accretion, and confirm the results 
from Broug h ct al. (2002), is to measure the structural pa- 
rameters (size and surface brightness) of BCGs with respect 
to the properties of their host cluster environment. 

Many authors have studied the eff ects of mergers on 
galaxy surface brightness profiles since iToomre fc Toonir3 
il972h suggested that dynamical friction, and the conse- 
quent merging, of pairs of disc galaxies could form elliptical 
galaxies. Numerical simulations indicate that the de Vau- 
couleurs r * surface brightness profile llde Van couleursll94Sh 
is rob ust to the effects of mergers (lBarneslll988l: INavarrol 
I1990I) . The profile of the remnant is only disturbed with re- 
spect to the de Vaucouleurs prediction at radii greater than 
the effective radius and these disturbances move outwards 
with time, disappearing within the dynamical timescale 
of the gala:xy. 'Schwcizcr' (Il982h observed that the V^-band 
light distribution of NGC7252, a galaxy which appears to 
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have undergone a recent merger, had already estabhshed 
an ri profile . IWrieht et~all (llOQCT l: iJames et al.l ^1999^ and 
iBrown et al.l (|200,?i ) examined A'-band surface brightness 
profiles of merger galaxies and merger remnants, all confirm- 
ing that the structure of the mergers more closely resembled 
that of elliptical galaxies than spiral bulges. 

More recently, mergers between elliptical galaxies have 
been simulated in order to understand the observed tilt 
of the Fundamental Plane (FP; between radius, surface 
brightness and velocity dispersion) with respect to the 
Virial theorem. These simulations find that the effective ra- 
dius of merge r remnants increases with total mass - (e.g . 
Navarrdll99C:ICapaccioU et alJll992l: ICapelato et al.lll995t 
D^^^^^ 2003r^vstim^eTO^^lTl2004) . with the sur- 
face brightness, consequently, decreasing. Equal-mass dis- 
sipationless merging of simulated elliptical galaxies pro- 
duces structurally non-homologous remnant s and is, there- 
fore, able to reproduce the observed FP iCaoelato et alJ 
ll995l:lNipoti et al.ll20o3) . These studies suggest that galax- 
ies become larger and more diffuse as a function of merger 
events. 

Observations of the Kor mendy relat i on (b etween ra- 
dius and surface brightness; iKormen dvl I1977I) of BCGs 
show that they lie off the relationship followed by nor- 
mal cluster elliptical galaxies, having larger radii and lower 
mean surface brightnesses than is predicted by simply scal- 
ing the relationshi ps followe d by normal elliptical galaxies 
iThuan fc Romanishin.l98U:lHoessel et aU1987: SchomberU 
ll987l:lQegerle fc Hoessellll99ll) . However theydo lie on the 
same FP as normal cluster ellipticals llOegerle fc Hoessel 
I1991I) . These observations are consistent with the suggestion 
that BCGs have undergone an increased merger history with 
respect to normal cluster elliptical galaxies. There are also 
suggestions that the de Vaucoule urs effective radii of BCGs 
correlate with their environment. ICarilli et all il997ll fitted 
de Vaucouleurs laws to the surface brightness profiles 
of 17 BCGs with redshifts z < 0.2 in the g—, r— and i— 
bands. They find a correlation between effective radius and 
local galaxy density (measured as the number of galaxies 
brighter than AIv = —18 in circles of projected radii equal 
to 200 kpc) in the sense that larger BCGs are found in en- 
vironments which are locally denser at the present epoch. 
Few other studies have taken environment into account so it 
is unclear what role environment plays in the structural dif- 
ferences observed between normal cluster elliptical galaxies 
and BCGs. 

In this paper the structural properties of BCGs are 
measured to further examine how the properties of BCGs 
depend on their host cluster environment, with the aim of 
interpreting their evolutionary history. This analysis uses a 
larger sample of objects covering a significantly wider range 
in environmental density than previous studies. An objec- 
tive measure of environmental density is vitally important 
in any study of BCG properties. As the cluster X-ray lumi- 
nosity is proportional to the volume integral of the square of 
the density of the intra-cluster medium it provides a quan- 
titative measure of the environmental density of BCGs with 
which to analyse possible environmental dependences. We 
also make use of near-infrared imaging with its associated 
benefits of small evolutionary corrections and low Galactic 
extinction 



The sample is outlined in Section|5|along with a descrip- 
tion of the reduction of the data. The analysis techniques are 
presented in Section |H] The results obtained are introduced 
in Section |1] and discussed in Section |K| and the conclusions 
drawn from this paper are presented in Section |S| 



2 DATA 
2.1 Sample 

Follow ing on from the photometric analysis in lBroueh et alJ 
i2002h the aim was to construct a sample of BCGs in X-ray 
selected clusters with a wide range of X-ray luminosities. 
Structural analysis requires high resolution images and to 
achieve this our analysis is primarily concerned with the 
structur al properties o f BCG s at redshifts z ^ 0.1. 

The iBurke et all J2000l) sample consists of 78 BCGs 
with redshifts 0.0 < z < 0.8 in clusters selected from Ein- 
stein Medium-Sensitivity Survey (EMSS; [Cioia & Luppiiiol 
I1994I) and the Sere ndipitous High-reds h ift Ar c hival ROSAT 
Clust er (SHARC; iBurke et alJ Il997l. l2003l : iRomer et all 
l2Q0Cf) X-ray cluster catalogues. These galaxies were observed 
in the if-band with the Infra-Red CAMera 3 (IRCAM3) 
on the United Kingdom Infra-Red Telescope (henceforth 
UKIRT; 1 pixel = 0.281" giving a field of view of 72") be- 
tween 1994 and 1998. Eleven of these BCGs are at redshifts 
2: < 0.1 (all from the EMSS catalogue) and, therefore, suit- 
able for use in this analysis. 

In contrast to the analysis in iBroueh et al.l (1200211 it 
was no t possible to use th e 2MASS extended source cat- 
alogue jjarrett et alj|20Q(J) for this analysis. 2MASS only 
provides one surface brightness measurement and the two 
images it provides are not adequate for a structural analy- 
sis: The postage stamp images are scaled to the size of the 
galaxy and are too small to analyse to a reasonable radius. 
The images from the actual scans, have a scale of 1" /pixel 
which is too large to provide a meaningful analysis. 

Therefore, to extend this analysis a further 26 BCGs 
with re dshifts 0.04 < z < 0.1 selected from the Lynam cat- 
alogue (lLvnam|[l999l) were observed over 3 nights (14-16 
January 2001) in the K„-band with the QUick Infra-Red 
Camera (QUIRC) on the University of Hawaii 2.2m tele- 
scope (henceforth UH; 1 pixel = 0.189" giving a field-of- 
view of 193"). These BCGs are in clusters selected from the 
Northern ROSAT Al l-Sky X-ray cluster survey (NORAS; 
iBohrineer et al J 1200(11 . The UKIRT and UH data were re- 
duced in the same manner and a brief overview of this pro- 
cess is given below. 

As there are two galaxies in common between the UH 
and UKIRT samples (those in the clusters MS0301-7 and 
MS0904-5) this gives a total sample of 35 BCGs, all with 
redshifts z < 0.1. 



2.2 Data Reduction 

The data reduction was performed using the 
IRAF/DIMSUM (Deep Infrared Mosaicing Software 
UM) package. 

The individual frames were masked for bad pixels, dark 
subtracted and divided by the exposure time. A flat field 
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image was created by median combination of the object im- 
ages, or separate sky exposures if these were available, and 
applied to the object frames. Masking of cosmic ray events 
was performed on the flattened images before they were mo- 
saiced together, which completed the processing of those 
objects with sky exposures. Otherwise the mosaic, which is 
substantially deeper than the individual exposures, was used 
to create an object mask, which was then applied to the in- 
dividual images before they were median-combined to form 
a flat. The flattened images were then processed as above to 
create the final image. 

Observations of stars from the UKIRT faint standards 
list JCasah fc HawardenI ITooj) were used to calibrate the 
UKIRT photometry onto the UKIRT system assuming an 
extinction of 0.088 mag airmass"^, the median value for K- 
band observations at Mauna Kea. 

The UH images were calibrated using observations of 
standard s tars selected from the U KIRT extended faint stan- 
dards list (|Hawarden et alJl200J) . At least 7 different stan- 
dards were observed at least 3 times each night. With little 
variation over the 3 nights, the observations were combined 
to give an atmospheric extinction coefficient of 0.118 mag 
airmass"^ and a zero-point of 23.305 with a standard devi- 
ation a = 0.029 mag. 



3 ANALYSIS 

One of the major problems in measuring the structural 
properties of any galaxy type is syst ematic errors in the 
measurement of the sky levels (e.g. iGraham et al.l llQQfit 
iGonzalez et al J 120011 ) as even the smallest excess flux can 
swamp the galaxy flux in the outer annuli, and oversubtrac- 
tion results in a truncation of the profile. This problem is 
increased when analysing BCGs, due to their larger size and 
the fact that, locally, BCGs are obser ved to have sig:nifica, nt 
amounts of flux in their outer wings iGraham et alj[l996l) . 

In t his analysis the sky levels were calculated using SEx- 
tractor llBertin fc Arnoiitall99(l ). In this program the sky 
level is estimated by constructing a grid and estimating the 
background within each square with a sigma-clipped mean. 
The grid is then median filtered to remove the effects of large 
bright objects. A grid size of 64 pixels and a median filter 
size of 3 X 3 grid squares was used here in order to avoid 
the background estimation being affected by the presence of 
other galaxies. This method was compared with measuring 
the background by obtaining the median of the counts mea- 
sured in small annuli at different points around the image, 
the difference was of the order of 10 per cent in (Jbackground 
calculated by SExtractor and no systematic offsets were ob- 
served. 

Owing to their position at the centres of clusters of 
galaxies, BCGs are often surrounded by many other galaxies. 
It is important not to incorporate the light of these galaxies 
when measuring the surface brightness profile of the BCG. 
In order to avoid this problem, the flux from the neighbour- 
ing galaxies needs to be carefully excluded from the images. 
When fitting surface brightness profiles in 2D it is possi- 
ble to deconvolve each galaxy and subtract the flux due to 
the curves of growth of the companion galaxies, replacing 
the actual BCG flux in each excluded pixel with that cal- 
culated using its symmetric counterpart about the centre of 



the BCG (e.g. lLaueilll98^ : iGarilh et alJll997l: iNelson et alJ 

■2002). However, this assumes that the BCG is symmetric 
and that there is not a further object affecting the counter- 
part position. 

To reduce systematic uncertainties in this analysis these 
assumptions were not made. Pixels that were identifled as 
belonging to a contaminating source were excluded from the 
profile calculation. The exclusion radius around each source 
was initially set to twice the semi-major axis of the contam- 
inating source (as reported by SExtractor). These apertures 
were then inspected by eye, and the exclusion radii of those 
sources close to the BCG centre were reduced, since the 
BCG flux close to the centre is so high that flux from other 
galaxies is insignificant. 

The galaxies were then fitted with elliptical isop hotes 
usin g the IRAF task ELLIPSE ('jedrzei ewskil Il987l: e.g. 
Gari lli et alj Il99^ : 'Gonzalez ot al. 2001) which calculates 
the median intensity within each isophote, further reduc- 
ing the effect fiux from companion galaxies can have on the 
profile of the BCG. The centres of the ellipses were fixed 
while the position angles and ellipticities were allowed to 
vary freely. 

ELLIPSE was allowed to run until either the angular 
extent of an elliptical annulus reached the width of the image 
or greater than 75 per cent of the pixels in a single annulus 
were masked out. Intensities differing from the median by 
more than 2 standard deviations were excluded to ensure 
the quality of the data being fitted. The radii were then 
circularized (r = y/ab) and converted from arcseconds t o 
h^^ kpc assuming a WMAP cosmology llSpergel et alj2do3l : 

= 0.3, D.A = 0.7 and Ho = WOh km s"^ Mpc"^ with 
h = 0.7. 

3.1 Petrosian Structural Parameters 

Petrosian structural parameters jPetrosianl Il976ll have a 
number of advantages over those calculated from fits of mod- 
els (e.g. de Vaucouleurs law, Sersic law) to the surface bright- 
ness profile. In particular, no assumption of the stellar light 
distribution is made in their calculation, making them a di- 
rect measure of galaxy properties, from the local surface 
brightness. The Petrosian parameters are also relatively in- 
sensitive to zero-point errors and errors in the extinction 
correction, as is clear from the Petrosian relationship: 

r|{r)=f^{^)-{^l{r)}, (1) 

where /i(r) is the surface brightness (in magnitudes) at a 
radius r and (/i(r)) is the mean surface brightness within r. 

Figure shows some example ri{r) profiles. The error 
bars indicate the la statistical errors on rj calculated by 
combining the measured errors on ^{r) and {^J.{r)) in quadra- 
ture. 

The Petrosian radii and surface brightness can be mea- 
sured from the rj profile at a chosen value of rj. The 
choice of ri value to use depends on the sample: r,, should 
be large enough to be una ffected by seeing (77 ^ 1.5; 
ISandaee fc Perelmuterlll9 90a') but not so large that the sur- 
face brightnesses are noisy and unduly sensitive to sky sub- 
traction errors. Figure [J illustr ates that the errors on the 
parameter increase with rj, and lKiaergaard et alJ (^9^ find 
larger errors in their at ry = 2 than at r; = 1.5. The images 
for this analysis have seeing of ~ 1.5" and the images have 
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Figure 1. Example ??(r) profiles of the BCGs in A1644, A644 (both high X-ray luminosity clusters), and those in MS0301-7 and A566 
(both low X-ray luminosity clusters). The caption on each panel gives the Petrosian radius, r,,, and surface brightness, /i,,, measured 
from each profile at r; = 1.5. The error bars indicate the Itr statistical errors on r) calculated by combining the measured errors on ^{r) 
and {fJ'ir)) in quadrature. 




log r (h ^ kpc) 



Figure 2. Example r)(r) profile for a model galaxy following a 
perfect de Vaucouleurs law {r^ = 5 kpc; fie = 17.8 mag 

arcsec"^). 



typical limiting depths of K~ 22 mag arcsec"^, therefore, 
Petrosian parameters calculated at = 1.5 satisfy both of 
these requirements. This results in a sample of BCGs with 
Petrosian radii in the range 7" < r^^i.s < 37" and sur- 
face brightnesses 19.5 mag arcscc^^ < fin=i.5 < 22.0 mag 
arcsec"^. 

The average Petrosian surface brightness (/i^) was cal- 
culated by summing the counts within and averaging 
over the area enclosed, (/x^) was then corrected for cos- 
mological dimmi ng (fl + zl**) and Ga lactic absorption us- 
ing the maps of ISchleeel et all (^9^; the absorption cor- 
rection is small, typically ~ 0.02 mag. The surface bright- 
nesses were also K-corrected using model predictions from 
the GISSEL96 co de jBruzual & Chariot 1993), following 
lYoshii fc Takaharal ^988). These K-corrections also include 
an evolutionary correction, assuming pure luminosity evolu- 
tion from a redshift of formation, zf — 2. The results are 
unchanged if the formation redshift is increased to 2/ = 5. 

The effects of masking companion galaxies and exclud- 
ing points with large errors in their intensities mean that 
it was not possible to calculate r,,=i,5 for 2 out of the 35 
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Figure 3. The distribution of cluster X-ray luminosity with red- 
shift for the host clusters of the 33 BCGs studied here. The 
squares indicate clusters selected from the NORAS X-ray cluster 
catalogue (Bohrineer et al. 2000), the triangles indicate those se- 
Iccted from the EMSS X-ray cluster catalogue iCioia &: Luppinol 
[l994). The dashed line indicates a cluster X-ray luminosity of 
Lx = 1.9 X lO^'^erg s~^, defined as marking the division between 
the two populations of BCGs at redshifts z > 0.1. 

BCGs in this sample, both in clusters with X-ray luminosi- 
ties < 1.9 X 10** erg s"'^. This leaves a total sample of 33 
BCGs with 0.04 < z < 0.1 in clusters with X-ray lumi- 
nosities 0.26 X 10*** < Lx erg s"^ < 5.8 x lO"*". The dis- 
tribution of the X-ray luminosities of the host clusters of 
these 33 BCGs is shown in Figure |21 This illustrates both 
the range of host cluster environment covered, over one or- 
der of magnitude in X-ray luminosity, and that this range 
remains consistent over the redshift range studied. The Pet- 
rosian structural parameters for these 33 BCGs are given in 
Table [U 

Statistical errors on the measurement of r,, and {^,,) 
were calculated by Monte Carlo realisations of the rj profile. 
The profile was resampled from a series of Gaussian profiles 
each with a FWHM equal to the rms of each point, and 
the parameters measured for each resampled profile. This 
process was repeated 1000 times resulting in 1000 indepen- 
dent measurements of each parameter, and, therefore, the 
spread around them. The standard deviation around the re- 
sampled parameters is of the order of 10 per cent in r,, and 
A{ij.ri) ~ 0.4 mag arcsec"^. 

3.2 Consistency 

It is possible to check the consistency of the analysis meth- 
ods described above as repeat observations were made with 
both telescopes. 

Table |21 summarises the Petrosian parameters mea- 
sured for the three BCGs with repeat observations made 
by UH and the two BCGs with repeat observations made 
by UKIRT. Two galaxies (those in the clusters MS0301- 
7 and MS0904-5) were observed by both UH and UKIRT, 
enabling a comparison of the photometry between the two 
telescopes. The 5 galaxies have a mean absolute difference 
Ar^ = 0.86 ± 0.27^"^ kpc and A(^^) = 0.20 ± 0.14 mag 
arcsec"^. The errors quoted are the error on the mean 



{— a/^/N). The difference in (/i,,) is comparable to the sta- 
tistical errors, however the statistical error on is signifi- 
cantly larger. 

These comparisons confirm that there are no significant 
offsets between the observations made on different nights or 
with the different telescopes. This also demonstrates the sta- 
bility of the parameters fitted and the validity of the analysis 
techniques used. 

4 RESULTS 

Two of the BCGs in Figure □ (those in Abell 644 and AbeU 
1644) clearly show a plate au in their r] profiles. This has 
been observed previously bv lKiargaard et alJlll993D in NGC 
4874, one of the two BCGs in Coma, which is known to 
have a cD morphology. This plateau is not present in the 
profiles of normal cluster elliptical g alaxies that follow per- 
fect d e Vauco ulcurs laws (Figure |21 S andaee fc Perelmuterl 
ll990allbl[T99ll:>Kiaergaard et al., 1993) . Ab ell 644 and Abell 
1644 h ave been classified as cD clusters bv lStruble fc RoodI 
il987f) . however, the BCGs themselves have not yet been in- 
dividually classified as cD galaxies. The plateaux we observe 
are analagous to those observed bv lKiargaard et al.l lll993fl 
in NGC 4874. We therefore speculate that the plateaux are 
likely to be a signature of the cD morph ology. The detailed 
study of the profiles of cD galaxies by ISchomberd (Il988h 
suggests that the envelopes of cD galaxies become domi- 
nant over the underlying galaxy component at an absolute 
surface brightness ~ ^ 24 mag arcsec"^, this is evi- 
dent as a break in the surface brightness profile at these 
surface brightnesses. Assuming a V — K colour ~ 3.0 mag 
from Bow er et al. (1992), this corresponds to hk ^21 mag 
arcsec"^ for these data. The plateaux in these data begin at 
r ~ 10 kpc and fix ~ 19 mag arcsec"'^ suggesting that 
the halo of BCGs with cD morphologies is distinguishable 
at higher surface brightnesses than found previously. 

The form of the Petrosian profile indicates a possi- 
ble fundamental difference between galaxies with cD mor- 
phologies and those that follow de Vaucouleurs profiles. The 
signal-to-noise ratio of these data preclude further analysis 
of BCGs with cD morphologies, and it is not studied further. 
However, it is noted that the Petrosian profile could provide 
an objective means with which to determine cD morphol- 
ogy without obtaining very deep imaging or fitting empiri- 
cal model s, which do not prov ide an adequate fit (Figure 
Figure 1^ ICraham et al.lll996ll . However, this requires fur- 
ther examination. 

The Petrosian Kormendy relationship is shown for all 33 
galaxies in Figure |1] The error bars show the statistical er- 
rors on these measurements. To illustrate any environmental 
dependence of BCG properties on the Kormendy relation, 
the cluster X-ray luminosity, Lx = 1.9 x l O^'^erg s~^ in 
the EMSS (0.3-3.5 keV) passband, defined in lBrough et alJ 
j2002i) as marking the division between the two populations 
of BCGs at redshifts z > 0.1, is used. There is an offset be- 
tween the properties of the BCGs on this relation, depending 
on their host cluster environment, but they both follow the 
same relationship. The best-fit Kormendy relation is given 
by: 

(m,) = 2.60'°-""' logio - 16.07(°'«^ . (2) 
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Table 1. The Petrosian structural properties. (^417) is corrected for K-dimming, cosmological dimming and Galactic absorption as 
described in the text. 
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0.27 


A1644 


0, 


,047 


2.073 


23.58 


2.21 


36.32 


19.63 


0.16 


A1775 


0, 


,071 


1.653 


29.27 


2.84 


30.61 


20.20 


0.26 


A1795 


0, 


,062 


5.797 


22.80 


3.42 


27.17 


20.06 


0.34 


A1800 


0, 


,075 


1.379 


18.94 


2.06 


18.89 


18.76 


0.19 


A1809 


0, 


,078 


0.940 


12.86 


1.71 


12.32 


18.65 


0.24 


MS0007-2 


0, 


,050 


0.487 


11.62 


1.50 


16.98 


18.58 


0.24 


MS0102-3 


0, 


,080 


1.154 


10.35 


1.55 


9.79 


18.11 


0.29 


MS0301-7 


0, 


,083 


0.293 


10.50 


1.78 


9.60 


18.91 


0.29 


MS0904-5 


0, 


,073 


0.829 


11.76 


1.63 


12.09 


18.32 


0.26 


MS1306-7 


0, 


,088 


1.494 


15.63 


4.01 


13.68 


18.64 


0.30 


MS1531-2 


0, 


,067 


0.405 


8.06 


1.00 


8.98 


18.44 


0.21 


MS1558-5 


0, 


,088 


1.249 


19.19 


1.62 


16.66 


19.38 


0.16 


MS1754-9 


0, 


,077 


1.315 


10.36 


1.75 


10.14 


17.97 


0.36 


MS2215-7 


0, 


,090 


1.049 


14.42 


2.27 


12.27 


19.32 


0.32 


MS2216-0 


0, 


,090 


1.697 


21.65 


2.34 


18.42 


19.39 


0.22 


MS2318-7 


0, 


,089 


2.688 


19.40 


1.66 


16.67 


18.95 


0.16 


MS2354-4 


0, 


,046 


0.372 


13.06 


0.85 


20.65 


18.41 


0.11 



Table 4. Statistical properties of the BCG Petrosian effective 
radii, r,, {h~^ kpc), and mean surface brightness, (/i,,) {K mag 
arcsec"^), as a function of cluster X-ray luminosity, Lx- 



Lx(xlO*Srgs-i) N r„ (la) (m^> (1<t) 

< 1.9 24 14.96 (6.23) 18.93 (0.68) 

> 1.9 9 20.03 (3.90) 19.33 (0.53) 



To examine the environmental dependence of these pa- 
rameters further the means and dispersions (split by the 
X-ray luminosity of their host cluster) are summarised in 
Table m 

Figure |S] illustrates the distributions of the BCG Pet- 
rosian radii (upper panel) and surface brightnesses (lower 
panel) as a function of their host cluster environment. A 
Spearman rank correlation on the Petrosian radii confirms 
that the BCG radii are correlated with their host cluster 
environment, with a correlation coefficient of = 0.55. For 
33 galaxies, the Student's t-test rejects the null hypothesis 
that there is no correlation at > 99 per cent confidence level. 
A weighted straight-line fit takes the statistical errors in r,, 



Table 3. Kormendy relations fitted to BCG Petrosian structural 
parameters, the relationships are of the form (fJ-rj) = A/xlogior,,-|- 
C. 



Author 


N 


A 


M 


C 


a 


Sandaee & Perelmuter fl991~) 


56 


V 


3.12 


8.51 


0.27 


JCiaergaar^^^al^ ^^93) 


18 


V 


2.11 


18.23 


0.31 


JCiajrgaar^^^al^ ^^22) 


18 


B 


2.09 


19.45 


0.29 


This sample 


33 


K 


2.60 


15.9 


0.42 


The quantities in brackets 


are 


the 


standard 


deviations 



around the gradient and intercept solutions. The scatter 
around the relation, a = 0.42 mag arcsec"^. 

The Kormendy relations fitted to other samples of Pet- 
rosian structural parameters are summarised in Table |21 
This sample follows a similar gra dient to those fit t ed by 
ISandaee fc Perelmutej Jl99ll) and iKiaergaard et alJ lll993lh 
to samples of BCGs, although it has a slightly larger scatter 
with a — 0.42 mag arcsec"^ in contrast to cr ~ 0.3 mag 
arcsec"^. 
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Table 2. Summary of the Petrosian parameters measured from repeat observations. The upper table shows the Petrosian radii, (h 
kpc), calculated for each night, the lower table shows the mean Petrosian surface brightness, (/irj) (mag arcsec"^). 



Cluster 


8 Nov 1994 9 Nov 1994 


{h 

11 Jan 1997 


-1 kpc) 
14 Jan 2001 


15 Jan 2001 


16 Jan 2001 


Mean 


A399 






25.99 


26.34 




26.17 


A978 






12.48 




12.76 


12.62 


MS0301-7 


11.51 








9.48 


10.50 


MS0904-5 


11.57 


12.47 




11.09 


11.89 


11.76 


MS1306-7 


15.74 


15.52 








15.63 



Mean Absolute Difference: 0.86±0.27 



Cluster 


8 Nov 1994 9 Nov 1994 


{^lr,) (mag 
11 Jan 1997 


arcsec ) 
14 Jan 2001 


15 Jan 2001 


16 Jan 2001 


Mean 


A399 






19.43 


19.43 




19.43 


A978 






19.02 




18.48 


18.75 


MS0301-7 


18.86 








19.09 


18.98 


MS0904-5 


18.35 


18.46 




18.29 


18.38 


18.40 


MS1306-7 


18.72 


18.69 








18.70 



Mean Absolute Difference: 0.20±0.14 



1^ . 

T- I ' ' ' ' — ' — r 




■^5 10 20 

(h"' kpc) 

Figure 4. The Petrosian Kormendy relation. The error bars 
indicate the Ict statistical errors on the measured parameters. 
The open points denote BCGs in clusters with Lx{0.3 — 3.5 
keV) < 1.9 X lO^^erg s'^ whilst the filled points denote BCGs in 
clusters with Lx(0.3 - 3.5 keV) > 1.9 X lO^'^'erg s-^ The solid 
line marks the best-fit relation described by Equation 1^ 

into consideration and also results in a statistically signifi- 
cant correlation between r,, and Lx'- 

rr, = S.TS'^-'^'logio Lx + 12.9<'-'>. (3) 

This is not a selection effect, Figure |21 illustrates that 
the host clusters of the BCGs have a uniform distribution 
in X-ray luminosity over the redshift range of the sample 
and Section 18.21 demonstrates that there are no offsets be- 
tween the parameters measured at the different telescopes. 
The BCG radii are clearly correlated with their host cluster 
environment. 

In contrast, the lower panel of Figure indicates that 
the Petrosian mean surface brightnesses are not strongly 



correlated with their environment. A Spearman rank corre- 
lation suggests that the Petrosian surface brightnesses are 
correlated with their environment with a correlation coeffi- 
cient Ts = 0.36. The Student's t-test rejects the null hypoth- 
esis of no correlation at > 95 and < 98 per cent confidence 
level, suggesting that there is a correlation of (fir,) with envi- 
ronment. A weighted straight-line fit demonstrates that the 
mean surface brightnesses of BCGs are weakly correlated 
with their host cluster environment. 

= 0.64(°"2) logio Lx + 19.0'°-^' . (4) 

The observed relationships between BCG parameters 
and those of their host cluster could be due to a system- 
atic error in the background subtraction of these galaxies. 
In order to test this, the Icrtackground calculated by SEx- 
tractor was used to model the possible systematic errors. 
This value is typically ~ 2 counts/pixel and is significantly 
more than the differences between the two background sub- 
traction methods compared in Section |3] The images were 
reanalysed - first with an extra IcJbackground subtracted and 
again with lobackground added. 

The Petrosian profiles of galaxies in those images with 
lo-background added fall off to zero before reaching rj = 1.5 
and it is not possible to measure their structural parameters. 
This is not seen in the original data and is illustrative that 
the systematic errors are not this large. 

The Petrosian radii of BCGs in images with 
lo-background Subtracted are smaller and the mean surface 
brightnesses are brighter than the original measurements. 
This is illustrated in Figure|S| which also demonstrates that, 
despite the significant background oversubtraction, the rela- 
tionships with the Petrosian structural parameters and the 
X-ray luminosities of their host cluster environments remain. 
The weighted straight-line fit to the upper panel of Figure|n] 

r, = 3.93(««'logioLx+8.l(0-«'. (5) 

The weighted straight- line fit to the lower panel of Figure |S| 
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O CM 



a. 

V 



L„ (10'*'* erg s * 0.3-3.5 keV) 



1 


1 

1 

,iJ 


n 1 1 1 — 1 — 1 — r-r-| 




1 


1 1 



1 



Lx (10'*'* erg s"* 0.3-3.5 keV) 



10 




L„ (10** erg s ^ 0.3-3.5 keV) 




Lx (10** erg s"' 0.3-3.5 keV) 



Figure 5. The distribution of Petrosian radius, (upper panel), 
and mean surface brightness {fir]) (lower panel), with cluster X- 
ray luminosity, Lx, for 33 BCGs. The error bars indicate the 
1(T statistical errors on the parameters. The dashed line indicates 
the empirical Lx = 1-9 x 10^'' erg s~^ division between high and 
low X-ray luminosity clusters. The solid line in the upper panel 
indicates the straight-line fit given by Equation El Equation 111 is 
illustrated by the solid line in the lower panel. 



(Mr,) =0.44(°-''''logio Lx + 17.6' 



,(0.2) 



(6) 



The BCG parameters are clearly correlated with those 
of their host cluster environment. 

The possibility that the observed relationships between 
the BCG parameters and those of their environment are 
due to a systematic error in the background subtraction 
of these galaxies w as further tested by a comparison with 
iGarilli et alJ (^oj) These authors examined BCGs in the r- 
band which is less susceptible to problems with background 
subtraction than the near-infrared, we have 7 BCGs in com- 



iGarilli et alJ il997ft fit de Vaucouleurs profiles to their 
BCGs and quote the fitted effective radius, r^, and effective 
surface brightness, He- Using these values, it is possible to 
construct the equivalent Petrosian profile for that de Vau- 
couleurs profile (as seen in Figure|5J. First, the effective ra- 
dius, re was transformed into the cosmology used here. The 
Petrosi an radius wa s then measured at rj = 1.5 as described 
above. iCarilli et alJ il997ll do not circularize their profiles 



Figure 6. The distribution of Petrosian parameters r,, (up- 
per panel), and mean surface brightness (/i,j) (lower panel) 
measured from images with the background oversubtracted by 
^^baakground- The error bars indicate the Icr statistical errors on 
the parameters. The dashed line indicates the empirical Lx = 
1.9 X 10*** erg s~^ division between high and low X-ray luminos- 
ity clusters. The solid line in each panel gives the straight-line 
fit fitted to the original data, and the dotted line indicates the 
straight-line fit to the oversubtracted data, described by Equa- 
tion |5] (upper panel) and Equation 151 (lower panel). 



as we have done here so these comparisons are made with 
the major-axis profiles of these galaxies. It was not, there- 
fore, possible to compare r,,(this work) with (comparison) 
for the BCG in MS1558-5 as the major-axis profile for this 
galaxy does not reach r] — 1.5. This leaves a sample of 6 
galaxies in common. The measured values are given in Ta- 
ble |S1 We also note that the of the BCG in Abell 1644 
measured by S. Andreon is in apparent agreement with this 
work (S. Andreon, private communication). 

The mean difference r7,(this work)— (comparison) = 
—9.3 ± 5.2 kpc. The most discrepant galaxy is the 

BCG in Abell 401, with r^fthis work)— r„ (c omparison) 
= 33.92 kpc. However, iGarilli et alj (^Q^) note that 
they determine a de Vaucouleurs effective radius for this 
BCG which is ~ 25 per cent larger than that obtained by 
iMalumuth fc Kirshneil l|l98,5h (re = 91.04 kpc). If we dis- 
count this galaxy the mean difference between the two sam- 
ples is then — 4.3±1.9 kpc. This difference is of the order 
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Table 5. Comparing major-axis r,, measured from de Vau- 
couleurs profiles fitted by Garilli et al. (1997) with galaxies in 
common with this work. 



Cluster 


re 

h"^ kpc 


Trj (comp) 

kpc 


(this work) 
kpc 


AVr, 

kpc 


A399 


35.12 


44.67 


36.20 


11.79 


A401 


59.67 


74.67 


40.75 


6.47 


A671 


21.74 


28.28 


23.92 


5.67 


MS0102-3 


9.41 


13.01 


13.97 


3.26 


MS0301-7 


10.99 


15.05 


13.92 


3.26 


MS0904-5 


15.66 


20.80 


12.24 


2.61 



of the errors on r^(this work) which are Ar^ ~ 4 kpc 
(Tablel^. We therefore conclude that the observed relation- 
ships are not due to a systematic error in the background 
calculation of these galaxies. 



5 DISCUSSION 

In order to understand what the observed correlation be- 
tween Petrosian radius and host cluster environment means 
in terms of BCG formation and evolution, the effect that 
mergers have on the surface brightness profiles of galaxies is 
considered. 

As introduced in Section equal-mass, dissipationless 
mergers of simulated elliptical galaxies produce structurally 
non- homologous remnants and are, therefor e, able to repro- 
duce the observed FP of elliptical galaxies iCapelato et alJ 
Il995l: iNipoti et al.l 120031) . However, these simulations are 
unable to reproduce the projections of the FP, the Faber- 
Jackson and Kormendy relations. The Faber- Jackson rela- 
tion is characterized by a lower velocity dispersion than 
predicted by the Faber- Jackson relation for a given lumi- 
nosity increase and the Kormendy relation by significantly 
large r radii than pred icted from observations I^ Navarrg 
Il990l: ICiotti fc van Albada 2001: Nipoti et al. 20Qa) . These 
simulations are, however, consistent with observations of 
BCGs. BCGs also lie on the FP for n ormal cluster ellip- 
tical galaxies JOeeerle fc HoesseJ IT99lh but are observed 
to have significantly lower velocity dispersions than pre- 
dicted by the Faber- Jackson relation for norma l clus- 
ter elliptical galaxies jMalumuth fc KirshneJ Il98ll . Il985t 
lOeeerle fc Hoessell Il99ll) and larger radii than predicted 
by the normal cluster elliptical galaxy Kormendy relation 
("Thuan fc Romanishin'lQSl; Ho essel et al. 1987: Sc hombertI 
fl987.:.Oegerle fc HoesseL199L) . This suggests that BCGs can 
form from equal-mass, dissipationless mergers of large lumi- 
nous elliptical galaxies. However, gas dissipation must be 
important in the formation of less massive elliptical galaxies 
in order to increase their central velocities, consistent with 
the observed Faber-Jackson and Kormendy relations. 

Numerical simulations of elliptical galaxy mergers in- 
dicate that the effective radii of merger remnants increase 
proportionally with total mass, by a factor o f 1 — 2 iHoesse 
198C : lNavarrol 1199 0: Capac cioh et alJll992l: Icapelato et al 
1995l:lNipoti et al.l2003i) . There is a range of increase factors 
because the properties of the remnant are dependent on the 
energy and angular momentum of the collision as well as the 



mass ratios of the progenitors. Head-on, equal-mass dissipa- 
tionless mergers produce the largest remnant s, with effective 
radii double those of their progenitors (e.g. ICapelato et alJ 
Qj)95; Nipoti et al. 2003). From the predictions of current 
simulations, the observation that BCGs in high X-ray lumi- 
nosity clusters are larger, by a factor of ~ 2, than those in 
low X-ray luminosity clusters at redshifts z; ^ 0.1 then sug- 
gests that these galaxies have undergone one or two major 
(equivalent mass) merger events, or several accretion events, 
more than those BCGs in \ow-Lx clusters. Figure^indicates 
that BCGs with cD morphologies have larger Petrosian radii 
than would be measured if they followed perfect de Vau- 
couleurs profiles. Therefore, the BCG radius - cluster X-ray 
luminosity relationship could also suggest that there are a 
higher proportion of BCGs with cD morphologies in high 
X-ray luminosity clusters than there are in low X-ray lumi- 
nosity clusters. Higher resolution images would be required 
to determine this in a non-subjective manner. 

Numerical simulations indicate that the density profiles 
of the remnant galaxies remain the same as those of their 
progenitors, therefore, the effective surface density should 
scale as le oc Mr/r^ |Navarro 1990). In magnitudes, this 
corresponds to a change in {^e) of up to 1 mag arcsec"^ in 2 
consecutive equal-mass mergers between elliptical gala xies, 
as predicted by the simul ations of lCapelato et alJ l|l993) and 
lEvstigneeva et al] i2004l) . The corresponding change in sur- 
face brightness predicted is within the observational scatter 
of this data, explaining the lack of a clear correlation of BCG 
mea n surface b rightness with host cluster environment. 

iBrough et a l. ( 2002) concluded that BCGs in high X- 
ray luminosity clusters showed little sign of an increase in 
mass since a ~ 1, in contrast to those in low X-ray lumi- 
nosity clusters which have increased in mass by up to a fac- 
tor of four over the same ti mescale. C ombining the results 
presented here with those of lBrough e t al. (2002) a consis- 
tent picture emerges. These two observations suggest that 
BCGs in high X-ray luminosity clusters assembled their stel- 
lar mass at redshifts > 1 and have been passively evolving 
since, in contrast to those BCGs in low X-ray luminosity 
clusters which are still in the process of assembling today. 
This is qualitatively consistent with the theories of hierar- 
chical structure formation which suggest that the less mas- 
sive systems (clusters and galaxies) we observe today as- 
sembled more recently than the most massive systems (e.g. 
■Kauffmann fc Charlot..l998.) . 



6 CONCLUSIONS 

In this paper, the Petrosian parameters measured from the 
surface brightness profiles of 33 BCGs in clusters covering 
a wide range in environmental density have been presented, 
with the aim of interpresting their evolutionary history. 

Analysis of the Petrosian r) profiles of this sample in- 
dicates that there is a deviation of the Petrosian profiles of 
some BCGs with cD morphologies from those of a galaxy 
following a perfect de Vaucouleurs law. This is suggested to 
be a signature of galaxies with cD morphologies and may 
indicate a fundamental structural difference. 

The structural properties of BCGs are clearly deter- 
mined by th e ir hos t cluster environment, as suggested by 
ICariUi et alJ (^93). BCGs in high X-ray luminosity clus- 
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ters are larger and fainter than those in low X-ray luminosity 
clusters. 

Research into other cluster galaxy populations suggests 
that their properties are driven, not by the global cluster 
environment, but by their local environment as measured by 
the density of their nearest neighbours, with those gala^xies 
in the least dense environments, showing the most signs of 
evolution. The unique position of BCGs means that their 
local environment is the global cluster environment and this 
affects their evolution such that those BCGs in the least 
dense environments have undergone more evolution since 
a ~ 1 than those in the most dense clusters. This provides 
an indication of how a single population of galaxies evolves 
hierarchically. It is likely that BCGs arise from the merging 
of satellites, initially drawn towards the cluster centre by 
coUimated infall along filaments in the early epochs of cluster 
formation, however, this work suggests that this process is 
continuing in the least dense clusters at the present time. 
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